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Abstract 

We consider a gravastar model made of anisotropic dark energy with an infinitely tliin spher- 
ical shell of a perfect fluid with the equation of state p = {1 — j)a with an external de Sitter- 
Schwarzschild region. It is found that in some cases the models represent the "bounded excursion" 



X 

' stable gravastars, where the thin shell is oscillating between two finite radii, while in other cases 

they collapse until the formation of black holes or naked singularities. An interesting result is that 
we can have black hole and stable gravastar formation even with an interior and a shell constituted 
of dark and repulsive dark energy, as also shown in previous work. Besides, in three cases we have 
a dynamical evolution to a black hole (for A = 0) or to a naked singularity (for A > 0). This is the 
first time in the literature that a naked singularity emerges from a gravastar model. 
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I. INTRODUCTION 

As alternatives to black holes, gravastars have received some attention recently par- 
tially due to the tight connection between the cosmological constant and a currently accel- 
erating universe J2| , although very strict observational constraints on the existence of such 
stars may exist j3|. ^ 

The pioneer model of gravastar was proposed by Mazur and Mottola (MM) [4] , consisting 
of five layers: an internal core < R < Ri , described by the de Sitter universe, an 
intermediate thin layer of stiff fluid Ri < R < R2 , an external region R > R2 , described 
by the Schwarzschild solution, and two infinitely thin shells, appearing, respectively, on 
the hypersurfaces R = Ri and i? = i?2- The intermediate layer is constructed in such 
way that Ri is inner than the de Sitter horizon, while R2 is outer than the Schwarzschild 
horizon, eliminating the apparent horizon. Configurations with a de Sitter interior have 
long history which we can find, for exarnple, in the work of Dymnikova and Galaktionov jsl- 
After this work, Visser and Wiltshire js] pointed out that there are two different types of 
stable gravastars which are stable gravastars and "bounded excursion" gravastars. In the 
spherically symmetric case, the motion of the surface of the gravastar can be written in the 
form 0], 

^R' + V{R) = 0, (1) 

where R denotes the radius of the star, and R = dR/dr, with r being the proper time of 
the surface. Depending on the properties of the potential V{R), the two kinds of gravastars 
are defined as follows. 

Stable gravastars: In this case, there must exist a radius Rq such that 

V{R^)=0, V'{Ro) = 0, V"{Ro)>0, (2) 

where a prime denotes the ordinary differentiation with respect to the indicated argument. 
If and only if there exists such a radius Rq for which the above conditions are satisfied, the 
model is said to be stable. Among other things, VW found that there are many equations 
of state for which the gravastar configurations are stable, while others are not |6|]. Carter 
studied the same problem and found new equations of state for which the gravastar is stable 
0], while De Benedictis et al [sl] and Chirenti and RezzoUa j^l investigated the stability 
of the original model of Mazur and Mottola against axial-perturbations, and found that 
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gravastars are stable to these perturbations, too. Chirenti and RezzoUa also showed that 
their quasi-normal modes differ from those of black holes with the same mass, and thus can 
be used to discern a gravastar from a black hole. 

"Bounded excursion" gravastars: As VW noticed, there is a less stringent notion of 
stability, the so-called "bounded excursion" models, in which there exist two radii Ri and 
i?2 such that 

\/(i?i) = o, y'(i?i)<o, \/(i?2) = o, i^'(i?2)>o, (3) 

with V{R) < for a G where R2 > i?i 

Lately, we studied both types of gravastars 
can indeed be constructed, although the region for the formation of them is very small in 
comparison to that of black holes. 

Based on the discussions about the gravastar picture some authors have proposed al- 
ternative models 15|. Among them, we can find a Chaplygin dark star 16|, a gravastar 



IG 



and found that, such configurations 



sure 



supported by non-linear electrodynamics [17|, a gravastar with continuous anisotropic pres- 



18| and recently, Dzhunushaliev et al. worked on spherically symmetric configurations 



of a phantom scalar field and they found something like a gravastar but it was unstable 19|. 



In addition, Lobo 



20| studied two models for a dark energy fluid. One of them describes 



a homogeneous energy density and the other describes an ad-hoc monotonically decreasing 
energy density, although both of them are with anisotropic pressure. In order to match an 
exterior Schwarzschild spacetime he introduced a thin shell between the interior and the 
exterior spacetimes. 



10 



ll| to the case where the equation 



In this paper, we generalize our previous works 
of state of the infinitely thin shell is given by p = {1 — 7)0" with 7 being a constant, the 
interior consists of a phantom energy fluid 20|, while the exterior is still the Schwarzschild 
space. We shall first construct three-layer dynamical models, and then show both types of 
gravastars and black holes exist for various situations. The rest of the paper is organized 
as follows: In Sec. II we present the metrics of the interior and exterior spacetimes, and 
write down the motion of the thin shell in the form of equation ([T]). In Sec. Ill we show 
the definitions of dark and phantom energy, for the interior fluid and for the shell. In Sec. 
IV we discuss the formation of black holes from standard or phantom energy. In Sec. V we 
analyze the formation of gravastar or normal star from standard or phantom energy. In Sec. 
VI we study special cases where we can not have the "bounded excursion". Finally, in Sec. 



VII we present our conclusions. 



II. DYNAMICAL THREE-LAYER PROTOTYPE GRAVASTARS 



The interior fluid is made of an anisotropic dark energy fluid with a metric given by 20 1 

dsl_ = -fidt^ + f^dr^ + r^d^^, (4) 
where dVL^ = dO"^ + sin^(6')ci(/)^ and 

= (1 + 6r^)^(l + 26r2)'^, 

where a; is a constant, and its physical meaning can be seen from the following equation (jH]). 
Since the mass is given by m(r) = 6r'^/[2(l + 26r^)] then we have that 6 > 0. The corre- 
sponding energy density p, radial and tangential pressures Pr and pt are given, respectively, 
by 

fujh\ ( 3 + 26r2 



\87T J V(l + 26r2)2^ 



vry V (1 + 2fer2)2 y 327r [(1 + 26r2)3(l + 6r2)] 
(1 + w)(3 + 26r2) [(1 + 3w) + 26r2(l + u)] 
-8a;(5 + 26r2)(l + 6r^)|. (6) 

The exterior spacetime is given by the Schwarzschild metric 

dsl = -fdv^ + f-^dr^ + r^dQ\ (7) 

where / = 1 — 2m/r — (r/Lg)^ and = a/S/A^- The metric of the hypersurface on the 
shell is given by 

dsl = -dr^ + R^{T)dn\ (8) 
Since ds'^ = ds"^ = we flnd that = = R, and 

fit' - f2R' = 1, (9) 

f^"-f = 1' (10) 



where the dot denotes the ordinary differentiation with respect to the proper time. On the 
other hand, the interior and exterior normal vectors to the thin shell are given by 



-R,i,0,0), 



< = i-R,v,0,0). 



Then, the interior and exterior extrinsic curvature are given by 

= ^(l + 6i?2)-"/2t|j4(l + 5i?2^-/25^2^2^2(l + 6i?2)-/2i^2. 



l + 26i?2 ' 
Kl. = K~, sin' (9), 



(11) 



{2bR'uj + 2bR' + 3u + 1) 



(1 + 2hR'Y^l + bRWP - (1 + 26i?2)"Vl + bRH'^ 
2(1 + bR'y/^l + 2bR')R'^ (1 + 2bR')-\l + bR^yHR + Rt - Rt, (12) 

(13) 

i{l + bR')R 



(14) 
(15) 



K^^ = v[{2Llmv + LjRR - L^Rv + R^v){2Llmv - LjRR - LIRv + R^i)) - 

2LlR'R?]{{2m - R)Ll + R^y\Llm - R^)L-^R-^ + Rv - Rv (16) 



KI, = -v{{2m-R)Ll + R^)Ll 



Since 



2l| 



[Kee] = K^oe ' K 
where M is the mass of the shell, we find that 



-M, 



M = v{2m - R) + 



i{l + bR')R 



l + 26i?2 • 

Then, substituting equations ^ and ( ITU]) into ( 12Ui) we get 

l + bR' + R'il + 2bR') 



R- 



1/2 



1 + 6/?2)-(c.+l)/4(l + 26i?2){.;+2)/2 • 



(17) 
(18) 

(19) 
(20) 



(21) 



In order to keep the ideas of MM as much as possible, we consider the thin shell as consisting 
of a fluid with the equation of state, p = (1 — 7)17, where a and p denote, respectively, the 
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surface energy density and pressure of the shell and 7 is a constant. Then, the equation of 



motion of the shell is given by 



2l| 



M + SnRRp = AnU'lT^pu'^n^] = AttR^ (T^f^W^n^, - Tl^uyf) , (22) 

where is the four- velocity. Since the interior fluid is made of an anisotropic fluid and the 
exterior is vacuum, we get 

M + SnRR{l - 7)cr = 0. (23) 
Recall that a = M/^AuR^), we find that equation (|23|1 has the solution 

M = kR^^''-^\ 



where k is an integration constant. Substituting equation 
rescaling m, b and R as, 

1 

m — )■ mk 27-3 ^ 
2 

b 6/c2^-3, 
R Rk~^, 

we find that it can be written in the form of equation ([1]), and 

V{R,m,Le,u,b,'y) = 



(24) 

into equation (!?!]) . and 



(25) 



— {bl^+^R^^'%^L 



h 2 _ /)2 

'2 "1 '^l 



- r^+X'^'^^'R^Ll - R^L%b\ + 2mL%R%\ + 
R%2b\ + /^^L^fof - 6^+2^2^-^"i?4^)/(6j/'i?4))^ + fef ~^"6^+ii?2L2 - 

^ 5-+25-|-+l^2^47-4 _ ^2a.+3^2^-.^2 ^ 5-+25-|-+l^2^2 ^ 

2bl^+\''RmLl - 26^+2^7^"+^mL2i? + bl^+X^R"" - 6^+26^^-+l^4) 

(26) 



where 



61 = 1 + 62 = 1 + 26/?2. 



(27) 



The exterior horizons are given by 



24| 



where y = (m/Lg)^, ip = arccos (3^/3y), r^h denotes the black hole horizon and denotes 
the cosmological horizon. Note that if y > 1/27 the quantity 3a/% is greater than 1, giving 
an imaginary angle ip. Thus, the horizons Vbh and Tc are imaginary and the spacetime 
becomes free of horizons. 

To guarantee that initially the spacetime does not have any kind of horizons, cosmological 
or event, we must restrict Rq to the range, 

Tbh < Rq <rh or Tc, (30) 

where Rq is the initial collapse radius. Clearly, for any given constants m, uj, b and 7, 
equation (l26l) uniquely determines the collapse of the prototype gravastar. Depending on 
the initial value Rq, the collapse can form either a black hole, a gravastar, a Minkowski, or 
a spacetime filled with phantom fluid. In the last case, the thin shell first collapses to a 
finite non-zero minimal radius and then expands to infinity. To guarantee that initially the 
spacetime does not have any kind of horizons, cosmological or event, we must restrict Rq to 
the range, 

Ro > 2m, (31) 

where Rq is the initial collapse radius. When m = = b, the thin shell disappears, and the 
whole spacetime is Minkowski. So, in the following we shall not consider this case. 

Since the potential (126|) is so complicated, it is too difficult to study it analytically. 
Instead, in the following we shall study it numerically. Before doing so, we shall show the 
classifications of matter, dark energy, and phantom energy for anisotropic fluids. 



III. CLASSIFICATIONS OF MATTER, DARK ENERGY, AND PHANTOM EN- 
ERGY FOR ANISOTROPIC FLUIDS 



Recently [23[, the classification of matter, dark and phantom energy for an anisotropic 
fluid was given in terms of the energy conditions. Such a classification is necessary for systems 
where anisotropy is important, and the pressure components may play very important roles 
and can have quite different contributions. In this paper, we will use this classification to 
study the collapse of the dynamical prototype gravastars, constructed in the last section. 
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TABLE I: This table summarizes the classification of the interior matter field, based on the energy 
conditions j^], where we assume that p > 0. 



Matter 


Condition 1 


Condition 2 


Condition 3 


Normal Matter 


p + Pr + 2pt>0 


p + Pr >0 


p + Pt>0 


Dark Energy 


P + Pr + ^Pt < (J 


P + Pr ^ U 


P + Pt >^ 






p + Pr <Q 


P + Pt>^ 


Repulsive Phantom Energy 


p + Pr + 2pt<0 


p + Pr >0 


p + Pt<0 






p + Pr <0 


p + Pt<0 


Attractive Phantom Energy 


p+Pr + 2pt>0 


P + Pr >0 


p + pt <o 






p + Pr <0 


P + Pt>0 



TABLE II: This table summarizes the classification of matter on the thin shell, based on the 

energy conditions j^]. The last column indicates the particular values of the parameter 7, where 
we assume that p >0. 

Matter Condition 1 Condition 2 7 

Normal Matter a + 2p>0 a + p>0 -1 or 

Dark Energy a + 2p < a + p>0 7/4 

Repulsive Phantom Energy a + 2p < a + p < 3 



In particular, we define dark energy as a fluid which violates the strong energy condition 
(SEC). From the Raychaudhuri equation, we can see that such defined dark energy always 
exerts divergent forces on time-like or null geodesies. On the other hand, we define phantom 
energy as a fluid that violates at least one of the null energy conditions (NEC's). We shall 
further distinguish phantom energy that satisfies the SEC from that which does not satisfy 
the SEC. We call the former attractive phantom energy, and the latter repulsive phantom 
energy. Such a classification is summarized in Table I. 

For the sake of completeness, in Table II we apply it to the matter field located on the 
thin shell, while in Table III we combine all the results of Tables I and II, and present all 
the possibilities. 

In order to consider the equations (jl]) and ([6]) for describing dark energy stars we must 



9 



analyze carefully the ranges of the parameter oj that in fact furnish the expected fluids. It 
can be shown that the condition p + pr > is violated for a; < — 1 and fulfilled for w > — 1, 
for any values of R and h. The conditions p + > and p + Pr + '^Pt > are satisfied for 
a; < — 1 and —1/3 < u < 0, for any values of R and b. For the other intervals of u the 



energy conditions depend on very complicated relations of R and b. See reference 



m. This 



provides an explicit example, in which the definition of dark energy must be dealed with 



great care. Another case was provided in a previous work [23|. Taking several values of u 
in the intervals — 1 < w < — 1/3 and > 0, we could not found any case where the interior 
dark energy exist. 

In order to fulfill the energy condition a + 2p > oi the shell and assuming that p = 
(1 — 7)0" we must have 7 < 3/2. On the other hand, in order to satisfy the condition 
0" + p > 0, we obtain 7 < 2. Hereinafter, we will use only some particular values of the 
parameter 7 which are analyzed in this work. See Table II. 

IV. STRUCTURES FORMED 

Here we can find many types of systems, depending on the combination of the constitution 
matter of the shell and core. Among them, there are formation of black holes, stable and 
"bou.ded exc.s,o,. ^^^^s, as ,t has a.eadv shown previous wo*s Q-Q, and 

even a naked singularity constituted exclusively of dark energy. All of them are listed in the 
table HI. 

As can be seen in the figures [H [3] and O depending on the value of the cosmological 
constant, we can see that V{R) = now can i have one, two or three real roots. Then, 
we have, say, Ri, where Ri+i > Ri. For = Li (corresponding to A = 0) If we choose 
Rq > R3 none structure is allowed in this region because the potential is greater than the 
zero. However, if we choose R2 < Rq < R3, the collapse will bounce back and forth between 
R = Ri and R = R2- This is exactly the so-called "bounded excursion" model mentioned 
in j^, and studied in some details in l^-[l^. Of course, in a realistic situation, the star 
will emit both gravitational waves and particles, and the potential will be self-adjusted to 
produce a minimum at R = R^tg Ur where V {R = Rstatic) = Q = V {R = Rstatic) whereby a 



gravastar is finally formed 
collapse of the shell. 



10l-|l2l|. For Rq < Ri a. black hole is formed in the end of the 
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The scenario above can significantly be changed if we consider A > 0. In this case for 
Le > Lc, we also have bounded excursion gravastars if R2 < Ro < Rs- However, for i?o < Ri 
the final structure can be now a black hole or a naked singularity since the presence of the 
cosmo logical constant above a certain limit (L*) eliminates the event horizon (its radius 
becomes imaginary), as can be seen in the tables IV, V and VI. This is the first evidence 
of a naked singularity formation from a gravastar model. Moreover for Le = Lc, then 
-R2 = -R3, a stable gravastar is formed if Rq = R2, while for Lg < Lc there is only one 
real root. Note that for any value of Le > L*, a naked singularity is formed for small 
initial radius of the shell. This is already present in the de Sitter-Schwarzschild solution 
2^ . since is the exterior cosmological constant which allows to relax the inevitability of the 
horizon formation (differently from the Schwarzschild solution). The news here are that we 
have found a source (the shell with an interior anisotropic dark energy fluid) which can be 
matched with the de Sitter-Schwarzschild vacuum spacetime and, depending on the values 
of the exterior cosmological constant and the total mass, can represent a naked singularity. 

Thus, solving equation (I2T1) for R{t) we can integrate R{r) and obtain R{r), which are 
shown in the figures Ej H] and [6] for the case G. 



V. CONCLUSIONS 



In this paper, we have studied the problem of the stability of gravastars by constructing 
dynamical three-layer models of VW [gI , which consists of an internal anisotropic dark energy 
fluid, a dynamical infinitely thin shell of perfect fluid with the equation of state p = (1 — 7)0", 
and an external de Sitter-Schwarzschild spacetime. 

We have shown explicitly that the final output can be a black hole, a "bounded excursion" 
stable gravastar depending on the total mass m of the system, the cosmological constant 
Le, the parameter u, the constant b, the parameter 7 and the initial position Rq of the 
dynamical shell. All these possibilities have non-zero measurements in the phase space of 
m, Le, b, u, 7 and Rq. All the results can be summarized in Table III. 

An interesting result that we can deduce from Table III is that we can have black hole 
and stable gravastar formation even with an interior and a shell constituted of dark and 
repulsive dark energy (cases H and I). Still more interesting is the case G, represented by 
figures [H [3] and where for small radius of the shell we have no formation of a black hole 
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TABLE III: This table summarizes all possible kind of energy of the interior fluid and of the shell 
and compares the formed structures in the two gravastar models (Ae = 0, Ag > 0). The letters 
SG, UG, BEG, BH, NS and N denote stable gravastar, unstable gravastar, bounded excursion 
gravastar, black hole, naked singularity and none, respectively. 



Case Interior Energy Shell Energy Figures Structures (Ag = 0) Structures (Ag > 0) 



A 


Standard 


Standard 


SG 


SG/UG/BEG 


B 


Standard 


Dark 


BH 


SG/UG/BEG 


C 


Standard 


Repulsive Phantom 


BH 


BH 


D 


Dark 


Standard 


N 


N 


E 


Dark 


Dark 


N 


N 


F 


Dark 


Repulsive Phantom 


N 


N 


G 


Repulsive Phantom 


Standard E El [5] 


SG/BEG 


BH/SG/UG/BEG/NS 


H 


Repulsive Phantom 


Dark 


BH 


BH/SG/UG/BEG 


I 


Repulsive Phantom 


Repulsive Phantom 


BH 


BH 


J 


Attractive Phantom 


Standard 


SG/BEG 


BH/SG/BEG 


K 


Attractive Phantom 


Dark 


BH 


BH 


L 


Attractive Phantom 


Repulsive Phantom 


BH 


BH 



(for A = 0) and a naked singularity (for A > 0). This is the first time in the literature that 
a naked singularity emerges from a gravastar model. Besides, the figures El H] and [6] give us 
examples of the dynamical evolution of a gravastar to a naked singularity. 

Finally, the opposite final fates, because of the absence or the presence of a positive 



25j that can exist a connection 



cosmological constant, reinforces the hypothesis proposed in 
between naked singularities and some kind of weakness of the gravitational field, compared 
to that associated to black holes. 
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FIG. 1: The potential V{R) and the energy conditions EC1= p + pr + 2pt, EC2= p + pr and 
EC3= p + pt, for 7 = -1, cj = 1.5, b = 0.01 and rric = 0.2097639045. Case G 

TABLE IV: This table show the calculated horizons using the equations [28] and [29l The symbol 
i = denotes the imaginary constant. See figure [TJ 



m 



Cosmological Horizon 



Black Hole Horizon 



0.2097639045 



0.8 

L*e 

1.2 

1.3104314 
1.5 
10 

oo 



0.4795669330+0.2235024144 i 

0.8577020327 
0.9974942529 

1.213250267 

9.783239233 

oo 



0.4795669330-0.2235024144 i 

0.5136246691 
0.4866382631 
0.4638970952 
0.4202701210 
2m=0. 4195278090 
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FIG. 2: These figures represent the dynamical evolution of the shell to a naked singularity for the 
potential given by the figure [U assuming Rq = R{0) = 0.7. 

TABLE V: This table show the calculated horizons using the equations [28] and [29l The symbol 
i = denotes the imaginary constant. See figure [3l 



m 



Cosmological Horizon 



Black Hole Horizon 



0.3775 



1.5 



LI 



1.97 
2.0053430150 
3 

20 

oo 



0.8943821768+0.3869863502 i 

1.1976796420 
1.2948106450 
2.5081539280 
19.61123828 
oo 



0.8943821767-0.3869863502 i 

1.0759952780 
1.0151227670 
0.8151917235 
0.7560805420 
2m=0.7550 



- Brazil (Ref. 2399/03). 
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FIG. 3: The potential V{R) and the energy conditions EC1= p + pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, w = 1.5, h = 0.01 and nic = 0.3775. Case G 

TABLE VI: This table show the calculated horizons using the equations [28] and [29l The symbol 
i = ^J—l denotes the imaginary constant. See figure [5j 



m 
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1.6 

12 
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oo 
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